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A model of AW UMa 
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ABSTRACT 

The contact binary AW UMa has an extreme mass ratio, with the more massive 
component (the current primary) close to the main sequence, while the low mass star 
at q « 0.1 (the current secondary) has a much larger radius than a main sequence star 
of a comparable mass. We propose that the current secondary has almost exhausted 
hydrogen in its center and is much more advanced in its evolution, as suggested by 
St§pieh. Presumably the current secondary lost most of its mass during its evolution 
with part of it transferred to the current primary. After losing a large fraction of its 
angular momentum, the binary may evolve into a system of FK Com type. 
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1 INTRODUCTION 

Since its discovery (Paczyriski 1964) , AW UMa has played a 
crucial role in our understanding of contact binaries by forc- 
ing all theories to explain both, its extremely low mass ratio 
of q ~ 0.08, and the apparent equality of effective tempera- 
tures of both components. Eclipsing binaries of W UMa type 
are in direct contact, sharing a common envelope around 
both components. Lucy (1968a,b) was the first to recog- 
nize that the convection of gas within a common envelope 
equalizes the entropy, and hence the effective temperature 
is approximately uniform and the colours of both compo- 
nents remain constant throughout the orbital phase. Yet, 
as the components have very different masses, most nuclear 
energy is generated in the more massive component (the cur- 
rent primary) , and it has to be redistributed throughout the 
common envelope. It was shown that this process is unstable 
on thermal (Kelvin-Helmholtz) time scale, which leads to re- 
laxation oscillations, with the matter being transferred from 
the more massive component to the less massive star, and 
vice versa (Lucy 1976, Flannery 1976, Robertson and Eggle- 
ton 1977, Yakut and Eggleton 2005, and references therein). 
If we assume, that the mass ratio reversal took place during 
previous evolution of such systems, the current secondary 
was the initial primary while the current primary was the 
initial secondary. Hereafter, we shall use "the primary" in- 
stead of "the current primary" and "the secondary" instead 
of "the current secondary". 

Hazlehurst (1970) was the first to suggest that contact 
binaries may be evolved off the zero age main sequence. 
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Stepieri (2003, 2006) suggested that the secondary is the 
more advanced in its evolution, in analogy with the "Algol 
paradox". It is not clear if this assertion is generally correct. 
We propose this is likely applicable to systems with extreme 
mass ratios, like AW UMa (Mochnacki and Doughty 1972), 
because the secondary's radius is so much larger than it 
would be on the main sequence. This "radius excess" is the 
primary reason to suspect that the secondary is very ad- 
vanced in its evolution. 



2 PARAMETERS OF AW UMA 

The following parameters of AW UMa were found on the 
Hipparcos website in Strasbourg: V = 6.90 mag, trigonomet- 
ric parallax tv = 15.13 mas, spectral type Fl. Assuming the 
bolometric correction BC = —0.1, these are combined to ob- 
tain the absolute magnitude of AW UMa: 

L = 6.61 L Q 

Assuming all the luminosity is due to the primary, our model 
on the Main Sequence gives 

Mi = 1.61 M s , X = 0.7, Z = 0.02, (by mass) 

where Mi denotes the mass of the primary component. This 
is consistent with a rather wide range of observational de- 
terminations of the total mass of the system (1.3 or 1.7 Mq 
(Rensing et al. 1985), 1.4 or 1.9 M Q (Rucihski 1992) which 
strongly depended on the assumed value of the mass ratio q 
in the interpretation of the radial velocity data. The binary 
period is 

Pbin = 0.4387 d 
For many years the mass ratio was adopted following 
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Mochnacki and Doughty (1972, Fig. 1) as q = 0.08. However, 
new high-quality data obtained at the David Dunlap Obser- 
vatory in 2006 (Rucinski, private communication) suggest 
that this value is too low and may be larger, q = 0.1 ± 0.02 1 . 
As the new result has not been published yet, we consider 
three values: q = 0.08, 0.10, 0.12. This choice of mass ratio 
slightly affects the size of the low mass secondary's Roche 
lobe, following Eggleton (1983). 



3 OUTLINE OF THE PROBLEM 

We propose a model in which the secondary was a star 
of ~ 1.5 Mq. It evolved off the main sequence, and it 
was stripped of most of its mass, down to the present 
0.14 — 0.18 M©. Some mass was transfered to the primary, 
some was lost from the binary. Also, some angular momen- 
tum was lost from the binary. We assume that a complicated 
evolution of AW UMa can be approximated with a model in 
which the primary has a structure of a single star, somewhat 
evolved off the Zero Age Main Sequence (ZAMS), while the 
secondary has a structure of a single star evolved up to the 
formation of its helium core, and stripped of most of its 
mass. In other words, we approximate the evolution of the 
two components of AW UMa with the evolution of two single 
stars. 



4 EVOLUTIONARY CALCULATIONS 

We adopt the initial chemical composition of AW UMa as 
X = 0.7, Z = 0.02, (by mass), and we use the evolutionary 
code as described in the readme file in: 

http://ftp.camk.edu.pl/camk/rs/04/readme.04 
This code follows evolution of a single, spherically symmet- 
ric star of a constant mass, in a standard manner. When 
necessary, the code was modified to take into account a 
rapid mass loss (Stage II, see below). As concerns input 
physics, we use the Livermore opacities (OPAL, Iglesias 
and Rogers, 1996) supplemented with molecular and grain 
opacities as given by Alexander and Ferguson (1994). We 
use, as well, the Livermore equation of state (Rogers et al., 
1996, Rogers, 2001). Nuclear reaction rates are calculated 
according to Bahcall and Pinsonneault (1995) updated ac- 
cording to Adelberger et al.(1998). This code has been al- 
ready used in some previous works (e.g., Dziembowski et al., 
2001). A grid of stellar models was calculated, with masses 
M 2 ,o/M© = 1.00, 1.28, 1.79, evolved from ZAMS until all 
hydrogen was burned out in their cores and helium cores 
were formed. This was referred to as Stage I in our evolution. 
During this standard evolution, the initial convective cores 
of the 1.28 and 1.79 M© models vanish and - during a short 
interval before the helium isothermal cores are fully formed 
and while in the radiative equilibrium - the rest of the hy- 
drogen is burned out. The initial mass of the secondary M2,o 
is the first parameter of our grid. We consider these models 
as possible progenitors of the secondary of AW UMa. The 
only constraint we put on the initial mass of the primary 

1 Dr. Rucinski writes that the DDO data indicate velocity 
field deviations from the contact model and will require special 
investigation. 



Table 1. Parameters of evolutionary advance in the initial models 
of the tracks shown in Figs. 1-3. X Ci o denotes central hydrogen 
content by mass. Mh c ,o denotes mass of a helium, i.e., mass of a 
hydrogen exhausted core. Values of all masses are in solar units. 
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Mi ; o is that it is significantly lower than M2,o to ensure that 
this star will be only slightly evolved when the secondary's 
radius reaches its Roche lobe. But we do not have to assume 
an extreme initial mass ratio of the components because a 
significant amount of initial mass of the secondary may be 
lost into the interstellar medium. 

A degree of exhaustion of hydrogen in the interiors of 
these three models, X Cj o, is the second parameter of our grid. 
We have only considered X Cj o as low as 0.0241 (by mass) or 
less (see Table 1). Adopted values of X c ,o are not explicitly 
shown in Figs. 1-3 to preserve their clarity. In these, for each 
line style representing a different M 2 ,o, the model tracks are 
shifted towards the left as the amount of hydrogen exhausted 
during Stage I increases. 

Next, for each M 2 ,o, we were stripping mass of the three 
evolved models assuming thermal equilibrium - the evolu- 
tion was frozen, i.e. there were no time dependent terms. 
This was Stage II, during which the mass was reduced to 
M 2 /M© = 0.18, 0.16, 0.14, consecutively. The value of M 2 
is the third parameter of our grid of models. This way, fol- 
lowing Stage II, we had 27 initial models, which were next 
evolved with no additional mass loss - this was Stage III. 
We looked at the radii of the models, expecting that some 
will be expanded enough to be acceptable as a model for the 
secondary of AW UMa. 

The 27 stellar model tracks are our guesses for what the 
secondary's evolution might be like. The primary star was 
just assumed to have the mass of Mi = 1.61 M©, and willing 
to accept matter from the secondary. The only possible way 
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Table 2. Exemplary sets of the initial "Stage II" parameters of 
AW UMa. Here qn denotes the mass ratio of the system and 
Pbin denotes the orbital period of the binary. See text for further 
explanations. 
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for it to influence evolution of the secondary is through the 
value of the secondary's Roche lobe radius, Rl2, being a 
function of Mi and M2, for a given binary period, Pbin- 

The mass and angular momentum loss from the system 
was, for simplicity of the argument, ignored. But we may 
interpret the whole process as analogous to the evolution of 
short-period Algols, as there exist among them objects with 
extreme mass ratios and very short orbital periods (e.g., R 
CMa: M 2 = 0.17 M ,Mi = 1.07 M ,P b in = 1.14 days; for 
more examples see Sarna and de Greve 1996, and Ibanoglu 
et al. 2006) . So the evolutionary scheme for AW UMa could 
be the one for a short-period Algol, such that a progenitor 
is a detached close binary which lost a fraction of mass and 
angular momentum via magnetized wind during its main 
sequence life so that its period is short enough at the begin- 
ning of Stage II. An additional mass and angular momentum 
loss takes place during nonconservative Stage II. Nonconser- 
vative mass exchange is accepted in many models of Algol 
formation (Sarna and de Greve 1996, Nelemans et al. 2000, 
Eggleton and Kiseleva-Eggleton 2002). 

In Table 2 we give a few examples of binaries in which 
the secondary is just filling its critical Roche lobe when hy- 
drogen in its core is nearly depleted (see Table 1). Adopt- 
ing reasonable values of the mass ratio from the range 0.3- 
0.8 we obtain values of the orbital period in the range 
Pbin = 0.89 — 1.37 days. This is just an example what the 
initial Stage II set of parameters might be. As one can see, 
to obtain the present value of the orbital period of AW UMa 
we have to assume further mass and angular momentum loss 
during the mass transfer process. 



5 DISCUSSION 

The most important result in this analysis are Figs. 1-3 
showing three possible scenarios of the secondary's compo- 
nent evolution during Stage III, after we have reduced its 
mass to the values M 2 /Mq = 0.18, 0.16, 0.14, respectively. 
Depending on the initial mass of the secondary M 2 ,o and 
its chemical composition at the end of Stage I, Stage III 
will start at various instants and continue for significantly 
different time intervals. 

We see that only for M 2 = 0.18 M Q (Fig. 1) does the 
secondary's radius R 2 exceed the secondary's Roche lobe ra- 
dius Rl2 (solid horizontal line) during the star's evolution, 
so the secondary is therefore able to overflow its Roche lobe 
and form a common envelope. This occurs for the highest 
value of M 2 = 0.18 M©. For masses M 2 0.16 M the sec- 
ondary will never overflow its Roche lobe (Fig. 2 and 3). 
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Figure 1. The time variations of the AW UMa secondary's 
radius for M 2 =0.18Mq. The different line styles represent 
its different initial masses of the secondary: solid lines for 
M 2 ,o = 1-79 Mq, dotted lines for M 2j0 = 1-28 M Q , and dashed 
lines for M 2j o = 1.00 Mq. For each line style three tracks are 
shown, corresponding (from right to left) to initial models more 
and more hydrogen exhausted during Stage I (see Table 1). Rl 2 
stands for the secondary's Roche lobe radius. 



This, combined with our estimation of the primary mass 
Mi = 1.61 M , requires the mass ratio for AW UMa to be 
q> 0.1. 

Such a star, having a very small mass for its radius is al- 
most an "empty" object. Its internal structure is represented 
in Figs. 4-6. We arbitrarily chose the initial mass of the sec- 
ondary to be M 2j o = 1-79 Mq and its present mass (i.e., the 
secondary mass) M 2 = 0.165 Mq, enough to ensure the for- 
mation of a common envelope. For a given binary period 
and assuming the primary mass Mi = 1.61 M the Roche 
lobe radius of the secondary is Rl 2 = 0.613 R . Solid lines 
in these figures represent internal structure of the model 
at the time when it has a maximum radius R 2 = 0.678 Rq 
during Stage III of evolution. Dotted lines show the internal 
structure after some time, in a more evolutionary advanced 
stage. 

When looking at the radius versus mass dependence 
(Fig. 4) we see that matter in the the secondary's interior is 
highly concentrated towards the center resulting in a very 
dense and small core and very extensive low density enve- 
lope. The effect is that the bulk of the secondary appears as 
if it was almost empty. 

The hydrogen inside the secondary is exhausted in its 
core (Fig. 5) and stays only in an outer shell. The thickness 
of the hydrogen layer will decrease during further evolution 
of the secondary. As we see in Fig. 6, nuclear energy gener- 
ation is limited to a very thin (in mass) hydrogen burning 
shell which becomes even thinner during further evolution 
of the secondary. 

Our models indicate that a model of AW UMa and simi- 
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Figure 2. The same as in Fig.l but for M 2 = 0.16 Mq. 
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Figure 3. The same as in Fig.l but for M2 = 0.14 Mq. 



lar to very low mass-ratio binaries involving a highly evolved 
secondary component encounters a limitation at the very low 
q end in that such secondaries cannot be made smaller than 
0.165 Mq or thereabouts. 

Further evolution of such systems will depend mainly 
on the evolution of the more massive component. The sec- 
ondary may become a degenerate dwarf and/or be consumed 
by the primary. In this scenario the system might evolve 
into a FK Com type star, as earlier proposed by Bopp and 
Rucihski (1981), Bopp and Stencel (1981). 
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Figure 4. An example of internal structure of the AW 
UMa secondary. In this figure, the internal mass-radius depen- 
dence is shown. The track parameters are, M2,o = 1-79 Mq, 
M 2 = 0.165 M , X c>0 = 0.014. The model with maximal radius 
during Stage III (denoted with "maximum") has been cho- 
sen to represent the internal structure of the low mass AW 
UMa component (the secondary). A slightly more evolved model 
("post-maximum") shows the direction of further evolutionary 
changes. For given binary period and assuming the primary mass 
Mi = 1.61 Mq and M 2 = 0.165 Mq the Roche lobe radius is 
R L2 = 0.613 R . 
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Figure 5. The same as in Fig. 4 but for the internal mass- 
hydrogen content relation. 
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Figure 6. The same as in Fig.4 but for the internal mass- 
luminosity relation. 
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